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Abstract 

The future deployment of cognitive radios (CRs) is critically dependent on the fact that the incumbent primary 
Q ' user (PU) system must remain as oblivious as possible to their presence. This in turn heavily relies on the fluctuations 

of the interfering CR signals. In this letter we compute the level crossing rates (LCRs) of the cumulative interference 
' created by the CRs. We derive analytical formulae for the LCRs in Rayleigh and Rician fast fading conditions. 

(N : 

, We approximate Rayleigh and Rician LCRs using fluctuation rates of gamma and scaled noncentral x processes 

CO ' respectively. The analytical results and the approximations used in their derivations are verified by Monte Carlo 

in 
o 

o\ 
o 



simulations and the analysis is applied to a particular CR aUocation strategy. 



Index Terms 

■ Cognitive radio, dynamic spectrum utilization, level crossing rates, average exceedance duration. 
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I. Introduction 

It is now well known [2], [3] that granting exclusive licences to service providers for particular frequency 
bands has resulted in severe under-utilization of the radio frequency (RF) spectrum. This has led to global 
interest in the concept of cognitive radios (CRs) or secondary users (SUs). These CRs are deemed to be 
intelligent agents capable of making opportunistic use of radio spectrum while simultaneously existing 
with the legacy primary users (PUs) without harming their operation. 

In addition to ensuring their own quality of service (QoS) operation, the most important and challenging 
task for the CRs is to avoid adverse interference to the incumbent PUs. Hence, it is necessary to develop 
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schemes that can help PUs avoid such harmful interference. The recently developed [1], [4] methods 
based on radio environment maps (REMs) [5], [6] can help achieve this goal very efficiently. In [4] only 
those CRs are allowed to operate in a particular time, frequency or space slot that do not reduce the 
PU signal to noise ratio (SNR) by more than some agreed penalty. However, the approaches in [1], [4] 
allow CRs to operate on the basis of average signal to interference plus noise (SINR) ratio values and 
do not consider the instantaneous temporal variation of the interference. Throughout the paper SNR or 
SINR values represent long term average values while the interference is considered on an instantaneous 
scale. Note that small scale variations in the composite CR interference can degrade the PU performance 
even though the CR levels may be acceptable on average. Thus, the determination of the rate at which 
the instantaneous interference crosses a particular threshold and the duration for which it stays above or 
below it, is an issue of core importance. For PU system designers, the following questions are important: 

• How large is the CR-PU interference and can it be controlled? 

• How often will the interference exceed a threshold? 

• How long does the interference stay above a given threshold? 

• How do these issues vary with the type of fading? 

These questions form the focus of this paper. In particular, we make the following contributions: 

• We determine the level crossing rate (LCR) and average exceedance duration (AED) of the CR-PU 
interference for Rayleigh and Rician fading channels, and various CR interferer profiles. 

• For Rayleigh channels, we approximate the LCRs using fluctuation rates of a gamma process. 
Similarly, for Rician fading we approximate the instantaneous aggregate interference with a fractional 
order noncentral variable to evaluate the LCRs. These approximations are validated via simulations. 

• For CR systems where the long term interference has an imposed maximum, results show that the 
LCR is maximum at or around the maximum interference threshold and is virtually zero 5 dB 
beyond this point. We also show that compared to Rayleigh fading, in line of sight (LOS) channels, 
the interference rarely crosses the threshold and when it does, it only exceeds the threshold value for 
a short duration. 

The rest of the paper is organized as follows: Section II characterizes the instantaneous interference and 
derives the LCR and AED results. In Section III we present simulation and analytical results. Finally, in 
Section IV we describe our conclusions. 



II. Instantaneous CR performance 



In any CR allocation policy, for example [4], even if the target SINR of the PU is exactly met the 
fast fading will result in fluctuations of the instantaneous SINR both above and below the target. As a 
first look at this problem we fix the PU signal power and consider the instantaneous variation of the 
interference only. Hence, in this section we focus on the instantaneous temporal behavior of the aggregate 
interference. For this purpose we evaluate the LCR (and thus the average exceedance duration (AED)) of 
the cumulative interference offered by the CRs. First we calculate the LCRs for a Rayleigh environment 
and then we characterize them for Rician fading conditions. 

A. LCRs for Rayleigh Fading 

For a given set of CR interferers, the instantaneous aggregate interference, lRay{t), is given by: 

N 

lRay{t) = Y,Ii\Ht)? (1) 

i=l 

where Jj represents the long term interference power of the ith CR, hi(t) is the corresponding normalized 
channel gain so that in Rayleigh fading \hi(t)\'^ is a standard exponential random variable with unit mean 
and N is the number of interfering CRs. Note that we fix the long term interference values, Ii, . . . , In and 
consider the variation of the fast fading terms, h^. From ([T]), the aggregate interference is represented as 
a weighted sum of exponential variables. Such weighted sums can be approximated by a gamma variable 
[7]. Simulated results show that the gamma fit is very good, but are not shown here for reasons of space. 
However, the corresponding LCR results are shown to be accurate in Figs. [il|3l Note that the exact LCR 
computation for such sums was given in [8] for the case of three and four branch maximal ratio combining 
(MRC) by providing special function integrals. Recently, more general expressions for arbitrary number 
of branches have been derived in [9]. However, the approach of [9] results in numerical difficulties, 
especially for large values of N, which can be the case for CR systems. Hence, an approximation is 
useful to overcome these problems and to provide a much simpler solution. Thus, approximate LCRs 
for (dl) can be found by calculating the LCR of the equivalent gamma process. The LCR for a gamma 
process has been calculated in [10]. Therefore, the crossing rate of lRay{t) across a threshold, T, can be 
approximated by: 




(2) 



where r = E{lRayit)f /Var{lRayit)), 9 = E{lRayit)) /Var^lRayif)) and i?(0) = pRay{Q) is the second 
derivative of the autocorrelation function (ACF) of lRay{t) at time lag, r = 0. Hence, to compute the 
LCR in Q only the mean, variance and ACF of the random process in ([T]) are required. The first two 
moments of O can be computed as E{lRay{t)) = J2iLi h and Var{lRay{t)) = J2iLi ^f- To calculate the 
ACF, note that: 

hi{t + t) = pi{T)hi{t) + ^(l~pi{T))ei it) , (3) 

where ei{t) is independent of hi{t) and statistically identical to hi{t). Assuming a Jakes' fading process, 
Pi{r) is the zeroth order Bessel function of the first kind, jQ{2Tr f^r) and is the Doppler frequency. 
Using ([3]) we have: 

N 

E[lR,y{t)lRay{t + r)] = ^ IJjE[\hi(t)\''\h,{t + r)n 

N , N 

= Y.^^h + [y^i}E\\hm\p\{r)Amf + (i - p\{r))\eAt)n 

N N N 

ij^j i=l 1=1 

^ N X 2 N 

^ j=l ^ i=l 

where in the second to last step above, we have used the fact that cross products have zero mean and that 
E[\hi{t)\^] = 2. The ACF of © is given by: 

_ E{lR^y{t)lR^y{t + r))-E{lRayit))E{lRayit + T)) 

' Vy<^r{lR^yit))Var{lRay{t + r)) 

and with the relevant substitutions, the ACF becomes: 



PRayVr) - AT ^ • W 



Finally, using the expansion Jo(2vr/£)r) = 1 — 7r^/|,r^ + . . ., the second derivative of the ACF needed to 
compute the LCR in Q is evaluated as: 



PRay{^) = -4vr^ ^=^ 7,^ - (V) 

Hence, the three parameters, r, 6 and -R(O), are available and Q gives the approximate LCR. 



B. LCRs for Rician Fading 



The instantaneous aggregate interference, Irndt), for this scenario is given as: 

N 

IRic{t)=Y,I^\h^{i)\\ (8) 

i=l 

where \hi{t)\ is Rician, with Rician K-factor denoted by A', and A^, Ji, I2, . . . , In are as defined in ([T]). 
Hence, Irndt) is a weighted sum of noncentral chi-square (x^) random variables. Note that standard LCR 
results for Ricians [11], [12] and noncentral variables [12] cannot be applied directly here. The work 
in [13] is for a single Rician and in [12] the LCR applies to the case where Ii = I2 = ■ ■ ■ = In and an 
exact noncentral arises with integer degrees of freedom (dof). Instead, using the same approximation 
philosophy as that used in the Rayleigh case, we propose approximating ([8]) by a single noncentral 
X^. This approach is less well documented but has appeared in the literature (see [14]). Also note that a 
scaled, rather than a standard, noncentral distribution is required for fitting and the resulting best-fitting 
distribution will almost certainly not have integer dof. A noncentral variable with v dof, non-centrality 
parameter A and scale parameter a has the following PDF: 

v-2 

= - exp (^^^ ^ J j (v^) , (9) 

where I{v-2)/2 is a modified Bessel function of the first kind with order [v — 2)/2. Fitting the PDF in 
dll) to the variable in ([8]) is performed using the method of moments technique so that the approximate 
noncentral x^ has the same first three moments as If(ic{t). The derivation details are outlined in Appendix 
I. Note that there can be numerical difficulties with the approach for certain values of h, ■ ■ ■ ,In- 
However, when this approach does not work it is straightforward to perform a numerical minimization 
of the difference between the true moments of the CR interference and the moments of the «— scaled 
noncentral x"^ variable. Values of A, v and a which minimize this difference can then be used. 

The LCR of a noncentral x"^ process with integer dof can be readily obtained from [12]. In particular, 
if we substitute R = T, a'^ = 1, M = v/2, = A and = fo in [12, Eq. (15)] we get the following 
expression for the LCR of the a— scaled noncentral variable 

LCR/^^^ = V^^fD{aT)^X^^e(^^h^{^/X^). (10) 

The result in (flOl) holds good for a noncentral process with integer dof. In Appendix II we show that 



this formula is also valid for non-integer dof. Note that a similar extension for a central with integer 
order [15] to a central with fractional order [10] has also been shown to be correct. 



C. AEDs 

We define the AED for both Rayleigh and Rician environments as the average time that the aggregate 
interference stays above a given threshold T [11]. Mathematically, 

where -F(T) gives the distribution function of the aggregate interference. Note that the exact CDFs of 
both lRay{t) and Irndt) can be found in [16]. 



III. Results 

In order to evaluate the accuracy of © and (flOl) it is important to use realistic values of Ji, /2, . . . , In- 
Hence, we use a particular CR access scheme [4] to provide these values. The decentralized selection 
algorithm in [4] employs a controller that considers CRs in their order of arrival. Each interferer is 
considered in turn and is accepted if the combined interference from previously accepted CRs and the 
current CR is less than some interference threshold. If a CR is not accepted, the next CR in the list 
is investigated. The Jj values are generated in [4] from randomly located CRs in a circular region and 
include path loss and shadowing effects. In [4] a threshold value is used which corresponds to the PU 
accepting a 2 dB loss in its SNR due to the presence of CRs. 

From 1000 simulations, using the above selection procedure two sets of interferers were selected. The 
first set selected had the highest variance. Only 3 CRs were accepted and there was a dominant interferer 
which accounted for 95% of the interference power. The second set had the lowest variance, representing 
the no dominant interferer scenario. Here, 18 CRs were accepted with the largest interferer only accounting 
for 16%. In addition to giving examples of engineering importance (presence or absence of a dominant 
interferer) these two cases also test the general applicability of Q and (flOl) over a wide range of interferer 
profiles. These sets were obtained using the following parameter values: shadow fading variance, a = 8.0 
dB, path loss exponent, 7 = 3.5, radius of PU coverage area, R = 1000 m, radius of CR coverage area, 
Rc = 100 m, CR density of 1000 CRs per square kilometer, an activity factor of 0.1 and fo = 25 Hz. 



LCR and AED of CR-PU Interference 

Figures [H [21 and [3] show the LCR (normalized by Doppler frequency) of the interference for different 
types of fading and interference profiles. The x-axis is also normalized by the rms value of the process 



interference (see Appendix I). Figure [T] shows the effect on LCR of increasing the Rician A'-factor, 
with the strong LOS case being considerably narrower than the non-LOS case. Figures |2] and |3] also 
show the value of the normalized interference threshold that restricts the long term average interference 
value in the CR system (as shown by the dotted lines). Note that there are multiple thresholds since the 
normalization is different for different channels. For all types of fading, the maximum LCR is observed 
close to this threshold value. This is because the CR allocation method gives a mean interference level 
close to the threshold. Even in strong LOS conditions {K = 10 dB), the interference shows a significant 
number of level crossings across the buffer due to the scattered component. Figure [2] shows the case of 
Rayleigh fading where the interference budget is dominated by a single large interferer with a number of 
smaller additional interferers. Also shown is the no dominant interferer case. Figure [3] shows the same 
results for a Rician channel with K = 10 dB. Figures |2] and |3] show that when there are many small 
interferers, the interference is more stable compared to the dominant interferer case. The results in Fig. |3] 
are quite promising. In near LOS conditions, the interference has a much lower level crossing rate across 
the interference buffer for the no dominant interferer case. Hence, it may be a desirable part of the CR 
allocation policy to avoid any single user which takes up a significant part of the buffer. Finally, for 
completeness, we show the AED results corresponding to Fig. [3] in Fig. HI As expected, the time spent 
by the interference above a threshold decreases as the threshold value increases. Therefore, for the no 
dominant interferer case, the interference seldom crosses the threshold (see Fig. |3]), and when it does, it 
only exceeds the threshold for a small period of time. Finally, we note that all figures show an excellent 
agreement between the analytical approximations and the simulations. 



In this paper we determine the LCR and AED for the CR-PU interference for Rayleigh and Rician 
channels. We have shown that LCRs in Rayleigh environment can be accurately approximated by LCRs 
of a gamma process. Similarly, while deriving LCR approximations in Rician conditions we have shown 
that the LCR of a noncentral process with non-integer dof has the same form as that of a noncentral 




where T is the interference power level and m2 is the mean-square 



IV. Conclusion 



process with integer dof. The LCR resuks show that it is desirable for the interference to be made up 
of several small interfering CRs rather than a dominant source of interference. The LCR of the former 
case is more stable than the latter. The AED results also show that the interference exceeds the threshold 
value for small periods of time in the latter case. 



Appendix I 

Let Y denote the random variable defined by (|9l). The first three moments of Y are [7]: 

E(Y) =a~\X + u) (12) 

E(Y^) = a-\{\ + iyf + 2{\ + iy) + 2\) (13) 

E{Y^) = a-%{X + uf + 6(A + z/)^ + 2A(A + z/) + 8(A + u) + 16A). (14) 

Similarly, suppose mi,m2 and m-^ denote the moments of Imcit) in ([8]) about origin. Expanding Imcit), 
Illicit) and Ifiicit) into multiple sums and taking expectation using standard results in [7] leads to: 

N 

mi = ^/, (15) 
1=1 



m2 

1=1 ' i=\ 

N \ 3 N N , / \ 2\ N 



1=1 ' i=l kjti,k=l ^ \ '/ /i=l^ \ < / \ < / / 

Now applying the method of moments, we solve = E{Y'') for k = 1,2,3 and obtain the following: 

A = 0.5a{am2 — ami ~ 2^i), y = ami — A, (18) 
where a is the solution to the following quadratic equation: 

a m^ — a m^ + Qam^ + 3a mim2 — 3a m^ + 8am2 — 8am-^ — 8mi = 0. (19) 



Appendix II 



In [12] a noncentral process, denoted r is considered. The only part of the derivation in [12] that 
requires integer dof is the proof that the conditional distribution of r given r is Gaussian with variance 
Var(r\r) = Aa^r where ct^ is a variance parameter. In this Appendix we show that this is true for a 
general noncentral process. The LCR of a stationary gamma process was first derived by Barakat [10] 
in an optics context building on previous results in [17]. This analysis is based on the representation [10] 



\E{x)\'^dx 



(20) 



A 



where A is the region of integration (an aperture in [10]), E(x) is a circular complex zero-mean Gaussian 
process and is the resulting gamma variable. If E(x) is allowed to have a constant non-zero mean then 
for certain A, the resulting Q has a scaled noncentral distribution with arbitrary degrees of freedom 
(not necessarily integer). For this noncentral case let E{x) = Ei{x) + jE2{x) and 



E^ix) + El{x)dx, 



(21) 



where Ei{x) and E2{x) are both non-zero mean Gaussian processes. The derivative of Vt is therefore 



tl= I [2Ei{x)Ei{x) +2E2{x)E2{x)]dx. 

J A 



(22) 



Now it is well known [8], [9], [12], [13], [18] that Ei{x), E2{x) are zero-mean Gaussian variables 
which are independent of Ei(x), E2{x) and each other. Let the distribution of both derivatives be denoted 
by A/'(0, cr^). Hence, conditioned on Ei{x) and E2{x) over x E A, the derivative, (l, is also zero mean 
Gaussian. The variance of (l conditioned on {Ei{x), E2{x)\x G A} is given by 



E 



E 



[ [ (2Ei(x)Ei(x) + 2E2ix)E2ix))i2E,{y)E,{y) + 2E2{y)E2{y))dxdy 

J A J A 

{AEl{x)El{x) + AEl{x)El{x))dx 



(23) 



since Ei{x) is independent of Ei{y) for x ^ y. Also, since E[Ei (x)] = a^, the conditional variance is 



4ct2 f [El{x) + El{x)]dx = 4a^n 

J A 



(24) 



Hence, (l has the representation (l = 2aQ^^'^Z where Z ~ A/^(0, 1) and (l has the conditional density 



Since (l ~ A/'(0, 4(T^f2), conditional on f2, the proof is complete. 
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Fig. 1. LCR results for different fading conditions with dominant interferers. The solid lines represent analytical results. Simulation values 
are shown by the circle, star and triangle symbols. 
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LCR results for the dominant and no dominant interferer cases in a Rayleigh fading scenario. The solid lines represent analytical 
Simulation values are shown by the circle and star symbols. The interference threshold values and their LCRs are shown by dotted 
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Fig. 3. LCR results for the dominant and no dominant interferer cases in a Rician (K = 10 dB) fading scenario. The solid lines represent 
analytical results. Simulation values are shown by the circle and star symbols. The interference threshold values and their LCRs are shown 
by dotted lines. 
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Fig. 4. AED results for the dominant and no dominant interferer cases in a Rician {K — 10 dB) fading scenario. The solid lines represent 
analytical results. Simulation values are shown by the circle and star symbols. 



